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Abstract. 

Since the discovery of its transiting nature, the super-Earth 55 Cnc e has become one of 
the most enthusiastically studied exoplanets, having been observed spcctroscopically and pho- 
tometrically, in the ultraviolet, optical and infrared regimes. To this rapidly growing data set, 
we contribute 42 days of new, nearly continuous MOST photometry of the 55 Cnc system. Our 
analysis of these observations together with the discovery photometry obtained in 2011 allows 
us to determine the planetary radius (1.990l55:^|o R®) and orbital period (0.7365417l^;^^2^^28 
days) of 55 Cnc e with unprecedented precision. We also followed up on the out-of-transit phase 
variation first observed in the 2011 photometry, and set an upper limit on the depth of the 
planet's secondary eclipse, leading to an upper limit on its geometric albedo of 0.6. 
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1. Introduction 

Currently, the most straightforward method for increasing the number of known tran- 
siting super-Earths around bright stars is to monitor the host stars of RV-discovered 
exoplanets with minimum masses in the super-Earth regime during their predicted times 
of inferior conjunction relative to the Earth. The planet 55 Cnc e is the first super-Earth 
whose transiting nature was discovered in this way. The host star of this system is a 
very bright {V = 5.95) G8 dwarf, and 55 Cnc c orbits it in only 0.74 days. The transit 
detection was achieved independently with the MOST space telescope (Winn et al. 2011) 
and with the NASA Spitzer infrared space telescope detection (Demory et al. 2011). 

Beyond the transit detection, both instruments have revealed further properties of 
this system and its innermost, transiting planet. The MOST light curve hinted at out-of- 
transit variations in phase with 55 Cnc e's orbit. The modulation, whose amplitude is too 
large to be due to reflected light from 55 Cnc e, may instead indicate the existence of tidal 
or magnetic interactions between the planet and its host star. Spitzer was employed to 
observe four occurrences of the planet's secondary eclipse using its 4.5 nm channel, from 
which a brightness temperature of 2360 ± 300 K was determined (Demory et al. 2012). 

I Based on data from the MOST satellite, a Canadian Space Agency mission operated by 
Microsatellite Systems Canada Inc. (MSCI; former Dynacon Inc.) and the Universities of Toronto 
and British Columbia, with the assistance of the University of Vienna. 
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The authors propose that the measured brightness temperature may suggest inefficient 
heat redistribution combined with a low albedo, but that it is also possible that the 
4.5 /im observations probe a deeper layer in the planet's atmosphere, thus potentially 
indicating the presence of a thermal inversion layer. Recently, Crossfield (2012) obtained 
NASA InfraRed Telescope Facility (IRTF) spectroscopic observations of 55 Cnc A and 
derived an improved flux density for the star, leading to a re-estimate of the Spitzer 
data-based brightness temperature of 55 Cnc e. From these measurements, he obtains 
a slightly lower temperature of 1950ti9o ^ ^"^^ planet, which expands the range of 
possible albedo values consistent with this temperature. 

In this proceedings, we describe the acquisition and reduction of new MOST obser- 
vations aimed at searching for the secondary eclipse of 55 Cnc e, verifying whether the 
previously observed phase variations persist, and improving the physical and orbital pa- 
rameters of this transiting super-Earth in Section [2j 



2. Observations and Data Reduction 

55 Cnc was observed nearly continuously using the MOST space telescope for about 
42 days from 2012 January 14 to February 25. The observations were acquired in Direct 
Imaging mode with an exposure time of 0.5 s per individual frame. The images were 
downloaded from the satellite in stacks of 80 frames, resulting in a total integration time 
of 40 s per data point. 

A raw light curve was extracted from the images using aperture photometry. The 
magnitudes were fitted to a linear function of sky background and x-position polynomials. 
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Figure 1. MOST 2012 light curve after de-correlation of magnitudes from sky background and 
X and y pixel position, and the running average correction for straylight variations at the MOST 
orbital period. 



The remaining routine correction is the filtering of straylight variations modulated with 



MOST Photometry of 55 Cancri 



3 



the satellite's orbit of 101.4 minutes. This was accomplished by dividing the light curve 
into 21 sections approximately 2 days long each, folding each section at the satellite's 
orbital period and computing a running average of this phased data. We then removed 
the resulting waveform from that section of light curve. The final light curve contains 62 
337 points and is shown in the bottom panel of Figure [T] 

3. Analysis of the Light Curve 

We used a Markov Chain Monte Carlo (MCMC) algorithm to analyze the MOST 
photometry with the goals of improving the values of the transit and orbital parameters 
of 55 Cnc e and searching for the out-of-transit variation originally observed in 2011 
(Winn et al. 2011) as well as the secondary eclipse. For model comparison, we use the 
odds ratio. The Bayesian Information Criterion (BIC) is used as an estimate of the 
marginal likelihood for each model (Carlin & Louis 2009). 

We fitted the following eight parameters: the orbital period (P), the mid-transit time 
(TO), the ratio of the planet to star radii (Rp/Ri,), the scaled semi-major axis [a/R-^), 
the orbital inclination (i), the depth of the secondary eclipse {5), the amplitude of the 
orbital phase modulation assuming a sinusoidal function with its minimum coinciding 
with the mid-transit time (a) and a normalisation factor corresponding to the flux just 
outside of transit (-Fnorm)- Including quadratic limb-darkening coefficients in our fit did 
not improve the odds ratio, suggesting the limb-darkening is not well constrained by our 
photometry. These two parameters were thus held fixed at ui = 0.648 and U2 = 0.117, 
the MOST bandpass values for the grid point of stellar parameters [T^ff = 5200 K, 
logg — 4.5, [Fe/H] — 0.2; A. Prsa, private communication) nearest to those of 55 Cnc. 
Gaussian priors were used for P and a/ R-i,, and uniform priors for the remaining fioating 
parameters. 

We carried out a MCMC analysis of the combined 2011 and 2012 light curves. The 
results of this analysis are listed in column 1 of Table [l] While we do not detect a 
significant non-zero secondary eclipse, we are able to use the best-fitted value of S and 
its uncertainties to place an upper limit on its depth. Because fitting for the secondary 
eclipse does not improve the model likelihood, we also ran the MCMC analysis omitting 
this parameter in order to obtain improved uncertainties on the remaining parameters. 
The results can be seen in column 2 of Tabled 

For all parameters but 6 the column of interest is therefore column 2. The combined 
photometry phased at the best-fitted period of 55 Cnc e is shown in the Figure [2] with 
the best fit model based on the parameters listed in column 2 superimposed. The orbital 
period we obtain is more precise than any published values for 55 Cnc e, and it is in good 
agreement with those values (Dawson & Fabrycky 2010). Our transit parameter values 
agree with those of Winn et al. (2011) within the uncertainties, though only marginally 
for i and a/R^,. In particular, our non edge-on best-fitted orbital inclination value has 
implications for the transit probabilities of the other planets in the system (Fischer et 
al. 2008), if their orbits are co-planar with that of 55 Cnc e to within a few degrees. 

4. Discussion 

4.1. The Radius and Composition of 55 Cnc e 

Our analysis of all the available MOST photometry of the system results in a planetary 
radius of 1.990tJ5:EJ^oi?e. Our new estimate, and the MOST 2011 (Winn et al. 2011) 
and Spitzer (Gillon et al. 2011) values agree with each other within their uncertainties. 
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Figure 2. 2011 and 2012 photometry phased at the orbital period of 55 Cnc e and averaged 
in 5-min phase bins. The red line is the best-fitting transit model without a secondary eclipse 
parameter, based on the values in column 2 of Table [T] See Section [3] for details. 



Table 1. Orbital and transit parameters for 55 Cnc e (2) 



Parameter 


2011 -1- 2012 photometry 
(with secondary eclipse) 


2011-1-2012 photometry 
(without secondary eclipse) 


Measured parameters 

Period - P (days) 
Mid-transit time - To " (BJD) 
Planet-to-star radii ratio - i?p/i?.* 
Transit depth - {B.p/R.t)^ (ppm) 
Scaled semi-major axis - a/B.^ 
Orbital inclination - i (deg) 
Amplitude of phase variation - a (ppm) 
Eclipse depth - 5 (ppm) 


2455962. 0693+oo!!i8 

070 + 29 

■^'^-30 
o t;24+0 041 
'^■°^*-0.042 

1 +1S 


n 7Q«t;41 7+0 0000025 

2455962.06971° 
01936+° 

^-29 
C-9O+0.042 

85.4+11 
34+12 


Derived parameters 

Transit duration - td (days) 
Semi-major axis - a (AU) 
Planetary radius - Rp (-R®) 


0.0662_o 0026 

0.01545_(, 00025 
1 08'^ + '' 079 
-L.300_Q Qg2 


n nccn+0 0035 
0.0660_o 0028 
n ni ej^c+O 00025 

1 qqn + ° 



We have determined the radius of 55 Cnc e (as measured in a broad optical passband) 
with a precision of 4%. We used the planetary mass derived by Demory et al. (2011) 
{Mp = 7.81 ± O.56M0) to obtain the mean density, and obtained a value of 5.489^Q gQ g 
cm~^, one of the best constrained for a super-Earth. The la upper limit on its density just 
reaches a composition in the coreless silicate regime (Elkins-Tanton & Seager), but the 
span of the uncertainty lies mainly within the volatile envelope planets area of theoretical 
mass-radius models. 

Gillon et al. (2011) examined two plausible compositions which match their mass and 
radius measurements for 55 Cnc e: an envelope consisting of 0.1% H/He by mass and 
one consisting of 20% H2O by mass. They argue, based on the calculations of Valencia et 
al. (2010), that the former would evaporate within a few million years. A water vapour 
envelope would also be depleted over time, but only very slowly; the process is expected 
to take several billion years. Given the estimated age of the star (10.2 ± 2.5 Gyr; von 
Braun et al. 2011), a water vapour atmosphere is by far the more likely scenario. Our 
value for the radius of 55 Cnc e is consistent with this model. 
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4.2. Brightness Modulation at the Planetary Orbital Period 

Winn et al. (2011) suggested the phase modulation discovered in the 2011 MOST pho- 
tometry, too large to be explained by scattered light from the planet, could be attributed 
to star-exoplanet magnetic interaction. Because the orbital motion of the planet and 
the stellar rotation are likely not synchronised (the best estimate of the stellar rotation 
period is ^43 days; Fischer et al. 2008) , we do not expect interaction manifesting itself 
as a region of star spots and chromospheric activity fixed in the co-rotating frame of 
the star at the sub-planetary point, but rather a travelling wave moving through the 
photosphere and chromosphere of the star. Such an effect can naturally be expected to 
have a changing amplitude, and in fact that is what Gillon et al. (2011) found in their 
re-analysis of the MOST 2011 photometry. 

When folding the entire 2012 MOST light curve, we do not find a phase variation at 
the period of 55 Cnc e as strong as that present in the 2011 light curve. This motivated 
us to inspect the 2012 data set more closely and, while always smaller than in the 2011 
data, the amplitude of the phase variation does appear to change over the course of the 
observations. Interestingly, if we fold in phase the photometry only between HJD 2455945 
(near the end of the portion of the light curve more affected by stray light modulation) 
and HJD 2455969, we find the amplitude of the phase variation (42 ± 16 ppm) is about 
twice what we find by folding the entire data set. This suggests that the amplitude of the 
variation varies on a timescale of weeks and perhaps even days, which is consistent with 
the finding of Gillon et al. (2011). This section of the light curve (phased and binned as 
before) is shown in Figure pj with the best-fitting model superimposed. 



1 .0004 _ 
1.0002 . 

X — 

^ 1.0000— ■■ 

> 0.9998 - 
■§ 0.9996^ 
0.9994 — 
0.9992 ^ 




-0.4 



-0.2 



0.0 

Orbital phase 



0.2 



0.4 



Figure 3. Approximately 2/3 of 2012 light curve phased at the orbital period of 55 Cnc e and 
averaged in 5-min phase bins. The red line is the best-fit transit model without a secondary 
eclipse parameter, and shows a mo re pronounced phase variation than when the complete 2012 
light curve is used (see Section [4!2| for details). 



4.3. Eclipse and Albedo Limits 
We have not detected a secondary eclipse for 55 Cnc e. Our best estimate for the depth 
of the eclipse is (5 = "1^22 PP™i ^-nd we can use this result to set an upper limit of ^ ^ 
17 ppm. 

We can express the ratio of the planet-to-star fluxes (Fp/F^) as a function of the 
geometric albedo Ag (the fraction of light reflected or scattered by the planet at optical 
wavelengths), planetary radius and semi-major axis: 



(4.1) 
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A geometric albedo of 1 (100% reflective) means the eclipse could be no deeper than 30 
ppm. Our upper limit on S leads to an upper limit {la) of 0.57 on the geometric albedo. 

Since wc do not have information on the planet's brightness variation throughout its 
orbit, we assume it has the reflective properties of a Lambertian surface. For a Lambert 
sphere, the Bond albedo is Ab = 1.5 * Ag. This relation leads to an estimate of 0.85 for 
the 1(7 upper limit on the Bond albedo of 55 Cnc e, in agreement with the results of 
Demory et al. (2012). 

If the surface of 55 Cnc e is molten and assuming the absence of an atmosphere (still 
possible if the planet is corclcss and made of pure rock) or the presence of a transmissive 
atmosphere, the planet's geometric albedo is predicted to be 0.6 (Kane et al. 2011). With 
our measurement of the geometric albedo, we can exclude this model with la confidence. 



5. Summary 

We have acquired 42 days of new MOST photometry of 55 Cnc. After a careful reduc- 
tion, we have used these data for several different types of studies. The new photometry 
has allowed the radius of 55 Cnc e {Rp — l.OOOi?®) to be determined to within 4%. 
Our result agrees with the hypothesis of Gillon ct al. 2011 that the planet likely has an 
envelope of high mean molecular weight volatiles such as water vapour and/or silicate 
clouds. Wc find that the phase variation observed in the 2011 MOST photometry persists 
but with an amplitude at most half of that measured in the discovery photometry. The 
data also allows us to place a limit on the depth of the secondary eclipse, and thus on the 
geometric albedo of 55 Cnc e of 0.57. This leads to an upper limit of 0.85 on the Bond 
albedo of the planet. With this value, we can exclude a highly refiective molten surface 
with la confidence. MOST photometry to be acquired in winter 2013 will bring this value 
into the realm of reflective silicate clouds, allowing us to place more constraining limits 
on the existence of such clouds in the atmosphere of 55 Cnc e. 
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